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ABSTRACT 


The  T-matrix  method,  a  technique  invented  by  Peter  Waterman,  (I <>65.  1969.  1971,  1977.  1980)  has  proven  to  be 
valuable  for  calculating  classical  scattering  from  bounded  targets  The  method  is  useful  for  acoustical, 
electromagnetic  and  elastic  scattering  and  is  based  on  coupled  boundary  integral  equations  *nd  partial  wave 
expansions  of  relevant  physical  quantities.  Upon  truncation  of  the  partial  wave  expansion,  which  are  in  principle 
infinite,  a  system  of  equations  k  derived  that  can  be  expressed  m  terms  of  finite  matrices  This  system  in  turn  leads 
to  an  expression  that  maps  the  incident  field  onto  the  scattered  field  The  outcome  is  a  powerful  algorithm  that 
yields  essentially  exact  results.  The  point  at  which  the  partial  wave  senes  is  truncated  is  always  an  issue  in  using 
the  method.  If  too  few  terms  are  included,  the  results  will  be  very  inaccurate  Conversely,  if  many  terms  beyond 
that  required  for  convergence  are  included,  at  best  the  calculation  will  be  needlessly  time  consuming  and,  ai  worst, 
can  lead  to  numerical  instability  due  to  numerical  round-off  errors.  What  would  be  ideal  is  a  technique  that  leads  to 
optimum  truncation,  that  is.  the  point  just  beyond  which  convergence  is  achieved  This  type  of  technique  can  be 
obtained  by  a  method  based  on  starting  with  a  “small*  T-macnx  and  in  which  we  build  up  the  dimensions  term  by 
term  via  simple  vector  operations  until  convergence  is  achieved.  The  algorithm  is  described  and  results  are 
illustrated  for  some  examples. 
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INTRODUCTION 

The  appropriate  number  of  basis  functions  must  be  used  when  generating  a  T-matrix;  if  too  few  are  used,  the  T- 
matrix  will  not  converge  (Waterman.  1965.  1969,  1971.  1977,  1980;  Wall.  1980,  Werby  and  Chin-Bmg,  1985; 
Werby,  1991).  When  generating  a  T-matrix  for  well-known  shapes,  a  few  rules  of  thumb  usually  pertain  to  the 
number  of  basis  functions  to  use,  but  not  for  arbitrary  shapes.  In  the  past  if  the  T-mainx  was  insufficiently  exact 
the  T-matrix  would  have  to  be  regenerated  with  more  (erms.  This  regeneration,  however,  did  not  guarantee  that  the 
new  T-matrix  would  converge.  With  this  in  mind,  an  iterative  algorithm  is  created  so  that  there  will  be  no 
guessing  as  to  how  many  basis  functions  is  needed  to  ensure  convergence.  This  paper  is  divided  into  the  following 
sections.  We  will  present  the  theory  behind  the  procedure  (starting  with  an  nth  order  T-matrix),  describe  what 
occurs  at  convergence,  and  show  some  numerical  results. 


THEORY 

The  iterative  procedure  is  based  on  the  method  of  bordering  (Faddeeva.  1959).  In  general,  we  start  with  a  (n  x  n) 
matrix  T.  Another  matrix  of  order  (n+l  x  n+1)  can  be  constructed  based  on  the  original  matrix  T.  It  is  thus 
possible  to  build  any  order  matrix  by  this  procedure,  since  the  expanded  matnx  T  order  (n+1  x  n+ 1 )  can  be  used  as 
the  next  starting  point  to  obtain  the  matrix  T  of  order  <n+2  x  n+2).  Let  us  then  begin  with  an  nth  order  T-matnx 
expressed  as. 
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with  a  Dirichlci  boundary  con&uon.  In  general,  if  we  expand  the  Tmatnx  of  order  in  x  ru  to  an  order  of  in-#  \  * 
n*l),  Eq.  (1)  becomes 
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where  r  is  a  row  vector  and  c  is  a  column  vector  and  %„+  j  -n4.i  is  a  scalar.  Now  Q‘ 1  is  expressed  as 
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where  P  is  an  (n  x  n)  matrix,  q  and  s  are  vectors  and  bn+l/i+i  is  a  scalar.  They  have  the  following  expression, 
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Finally  we  can  rewrite  Eq  (5)  as 
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Using  this  algorithm  has  several  advantages  As  already  mentioned  the  algorithm  takes  the 

guesswork  out  of  choosing  the  number  of  basis  functions.  The  T-matru  ts  conditioned  as  rt  is  built,  this  technique 

relaxes  some  of  the  complications  associated  with  ill -conditioning,  which  ts  sometimes  encountered  for  various 

objects.  Finally,  inspection  of  the  theory  shows  that  the  addition  of  additional  columns  and  rows  ts  an  operation 

involving  either  a  matrix  and  a  vector  or  two  vectors,  which  if  performed  on  a  vector  machine  would  be  fast  and 

economical. 


DISCUSSION  OF  CONVERGENCE 

Partial  wave  series  solutions  of  scattering  and  rearrangement  collisions  have  long  been  used  in  quantum  scattering 
There  are  two  main  reasons  for  doing  so.  First,  it  is  much  easier  to  solve  for  each  of  the  partial  wave  solutions, 
which  reduce  from  a  three-dimensional  problem  to  that  of  a  senes  of  one -dimensional  problems  for  central  forces,  or 
in  general  for  problems  in  the  which  angular  contribution  can  be  separated  from  the  radial  pan.  Second,  convergence 
is  generally  expected  at  some  value  of  angular  momentum  due  to  the  centrifugal  repulsion  effect  What  this  means 
is,  that  for  some  value  of  ka  where  k  is  the  wave  number  and  a  is  some  characteristic  dimension  of  the  object  (the 
radius  for  spheres,  for  example),  at  some  point  the  centrifugal  contnbubon,  which  occurs  in  a  partial  wave  solution 
(namely  l(l+l)/rz).  is  greater  than  the  component  (Va)^.  Thus  the  contribution  from  that  centrifugal  term  is 
repulsive  and  causes  little  overlap  with  the  target  when  it  dominates. 

Figure  I  illustrates  this  point  for  several  incident  partial  waves.  As  evident  from  the  figure,  when  the  partial  wave 
number  increases,  there  is  a  resulting  decrease  in  the  overlap  with  the  target.  For  suitably  low  values  of  I.  the  ka 
component  dominates,  so  the  interaction  region  is  significant.  The  question  for  what  value  of  I  the  centrifugal 
component  will  dominate  has  been  answered  for  spheres  by  using  the  WKB  approximation.  It  is  generally  agreed 
that  a  value  of  1=2  ka  is  large  enough  for  the  series  to  cut  off.  For  spheroids  the  situation  is  more  complicated, 
since  the  system  of  partial  wave  solutions  are  coupled  and  since  the  overlap  of  the  incident  partial  waves  with  the 
surface  displacement  terms  is  rather  complicated. 

For  the  spherical  case  there  is  mode  coupling  of  the  surface  boundary  terms  with  each  partial  wave,  and  this  corrupts 
the  simple  picture  we  have  for  spheres.  Figure  2  illustrates  the  case  of  convergence  for  several  spheroids  with 
varying  aspect  ratios  that  deviate  between  2  and  18  for  a  kL/2  fixed  at  6,  where  L  is  the  object  length  and  k  the  total 
wavenumber  (k=2nA).  We  use  a  low  accuracy  T-matrix  code  so  that  numerical  difficulties  will  be  accentuated.  For 
low  aspect  ratio  targets,  convergence  is  maintained  with  increasing  partial  wave  terms.  As  the  number  of  partial 
wave  terms  increase,  more  terms  are  required  for  convergence  because  of  an  overlap  effect  with  the  surface  expansion 
terms  (Werfcy  and  Chin-Btng,  1983;  Werby,  1991).  Further,  the  point  at  which  the  answer  is  no  longer  accurate  is 
due  to  ill-conditioning  that  results  from  high  aspect  ratio  targets  which  occurs  earlier;  thus,  the  "window  of 
convergence"  narrows  and  makes  the  problem  of  convergence  more  urgent  with  increasing  aspect  ratio.  The  method 
outlined  in  the  previous  section  not  only  obviates  the  problem  of  convergence  but  also  appears  to  condition  the 
matrix  as  well.  Alternative  methods  exist  that  allow  problems  to  be  solved  based  on  the  coupled  equations  of 
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Waterman  for  elastic  targets  (Werby  and  Green.  1983/  and  I  n>m  the  paint  of  vie*  ot  an  eigen,  j^r  <hv>.  ie:h  .  V.  eib. 
<t  al  .1986)  The  problem  of  a  target  in  a  wave  guide  has  been  discussed  extensively  in  a  recent  anule  b<  V-rur 
and  Werby  (199!) 
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Fig  !.  Comparison  of  the  overlap  with  the  target  for  several  incident 
partial  waves. 
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Fig.  2.  Form  function  vs.  matrix  order  for  various  aspect  ratio  targets. 
L  is  the  length  and  D  is  the  diameter  of  the  target  The  region 
between  arrows  represent  the  "window  of  convergence." 


APPLICATION  TO  RIGID  TARGETS 

There  are  two  classes  of  targets  for  impenetrable  problems,  i.e.,  soft  and  hard  scauerers.  They  do  not  support  body 
resonances;  therefore,  we  examine  acoustic  quantities  appropriate  for  such  non-resonant  targets  as  circumferentially 
diffracted,  or  creeping,  waves.  These  waves  arise  when  scattering  end-on  from  a  spheroid  in  which  one  the  return 
signal  is  observed  at  the  origin  of  the  signal.  There  are  two  competing  mechanisms.  One  arises  from  specular 
scattering  (geometric)  and  the  other  arises  from  the  creeping  waves.  The  result  is  a  coherent  effect  in  which  the  two 
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Fig.  4.  a)  Bistatic  scattering  from  30  to  1  aspect  ratio  rigid  spheroid 
end  on;  b)  broadside  incidence  for  kL/2=200. 


SUMMARY 

We  have  described  an  iterative  method  that  produces  an  (n+m  x  n+m)  main*  from  an  (n  x  n)  matnx  With  the  use 
of  an  illustration  we  described  what  must  occur  in  order  to  produce  convergence.  Finally,  we  presented  examples  of 
scattering  results  achieved  using  this  technique. 
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